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We employ the variational method for the Sturm-Liouville eigenvalue problem to analytically 
study the phase transition between the holographic insulator and superconductor in the Gauss- 
Bonnet gravity. By investigating the s-wave and p-wave holographic insulator/superconductor 
O ' models, we find that this analytic method is more efTective to obtain the analytic results on the 

^ ' condensation and the critical phenomena in the AdS soliton background in Gauss-Bonnet gravity. 

W)| Our analytic result can be used to back up the numerical computations in the AdS soliton with 

Gauss-Bonnet correction. 
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I. INTRODUCTION 



The profound finding of the anti-de Sitter/conformal field theory (AdS/CFT) correspondence l|-|3| has 
provided a framework to describe the strongly coupled field theories in a weakly coupled gravitational system. 
Recently, this correspondence has been applied to study the holographic model of superconductors in which a 
remarkable connection has been observed between the gravitational physics and the condensed matter physics 
J|. It has been shown that the bulk AdS black hole becomes unstable and scalar hair condenses below a 
critical temperature. The instability of the bulk black hole corresponds to a second order phase transition 
from normal state to superconducting state which brings the spontaneous U(l) symmetry breaking. In the 
boundary dual CFT, these properties exhibit the behavior of the superconductor 5|, |6|. Due to the potential 
applications to the condensed matter physics, the condensation in bulk AdS black holes has been investigated 
extensively, for reviews, see Refs. 7|-|9[ and references therein. 

In additional to the bulk AdS black hole spacetime, recently it was found that a holographic model can be 
constructed in the bulk AdS soliton background to describe the insulator and superconductor phase transition 
lo| . Adding the chemical potential to the AdS soliton, a second order phase transition can happen when the 
chemical potential is over a critical value. This phase transition can be used to describe the transition between 
the insulator and superconductor, while it is different from the Hawking-Page phase transition between the 
Ricci flat AdS black hole and the AdS soliton [11]. Taking the backreaction of the matter fields into account, 
it was argued that the order of the phase transition can be changed from the second to the first if the 

n 

backreaction is strong enough jl2:]. In the Stiickelberg mechanism, rich physics on the phase transition 



between the holographic insulator and superconductor in the AdS soliton background has been observed 



13|. 



Further investigations on various insulator and superconductor phase transitions in different theories of gravity 



have been carried out 
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18i. 



In most cases, the holographic superconductors were studied numerically. Ideally, one would like to have a 
full analytic description of the phase transition and condensation phenomena. In addition to back up numerical 
results, the analytic description can help to gain more insights, for example it may tell what properties of the 
action decide the mean-field behaviors etc. Recently, there appeared two analytic approaches in parallel to the 



numerical calculation. One is the analytic matching method which was first proposed in 19| and later refined 



20j. With this method, we can calculate the critical temperature analytically within a few percent in 



the best case. This analytic approach has been extended to derive the upper critical magnetic field when the 
holographic superconductor is immersed in constant external magnetic fields 



21 



22|. In higher dimensions, 

3- 



this analytic method can keep valid only when the matching point is chosen within an appropriate range 
However the matching method is not effective to describe the AdS soliton, neither can it be used to derive 
the critical exponents for the condensation near the critical temperature. The mean-field critical exponent 
1/2 at the critical temperature comes mostly from numerically solving the holographic systems and doing 
data fitting. In 



23 



24| . the authors extended the variational method for the Sturm-Liouville (S-L) eigenvalue 
problem to analytically calculate the critical exponent near the critical temperature. This method was further 
applied to analytically study some properties of holographic superconductors in AdS black hole backgrounds 
in Einstein gravity [25] and Einstein-Gauss-Bonnet gravity [26] in the probe limit, respectively. 

An analytic study by using the S-L method on the phase transition between the holographic insulator 
and superconductor was done in 17[ in the Einstein gravity. It is of interest to further generalize the S-L 
method to study holographic superconductor developed in the AdS soliton background in the Gauss-Bonnet 
gravity. The condensation phenomena and the phase transition between the s-wave holographic insulator and 
superconductor in the Gauss-Bonnet gravity were investigated numerically in [l^ . We will also extend 
the investigation of the p-wave holographic insulator and superconductor phase transition with Gauss-Bonnet 
correction in this work, which has not been constructed as far as we know. It is not trivial to analytically 
study the condensation and the phase transition by taking into account of the influence of the Gauss-Bonnet 
coupling. Besides to be used to check numerical computation, the analytic investigation can clearly disclose 
the critical exponent of the system at the critical temperature and the influence of the Gauss-Bonnet factor 
in the phase transition. In the AdS soliton background, we will compare two available analytic methods and 
argue that the S-L method is more effective for the analytic study of the condensation. 

The plan of the work is the following. In Sec. II we briefly review the AdS soliton background in the 
Gauss-Bonnet gravity. In Sec. Ill we explore the s-wave insulator and superconductor phase transition with 
Gauss-Bonnet correction. In Sec. IV we discuss the p-wave case. We conclude in the last section with our 
main results. 
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II. GAUSS-BONNET ADS SOLITON 

In order to study the superconducting phase dual to the Guass-Bonnet AdS sohton configuration in the 
probe limit, wc start with the five-dimensional AdS soliton in the Gauss-Bonnet gravity in the form 

dr2 



ds' = -r'dt' + TFT + /(^) V + r'idx' + dy'), 



(1) 



with 



(2) 



where rs is the tip of the soliton which is a conical singularity in this solution, a is the Gauss-Bonnet coupling 
constant and L is the AdS radius. It should be noted that in the asymptotic region (r — oo), we find 



(3) 



so the effective asymptotic AdS scale can be defined by 

,0 2a 



^cff 



L2 

1 



for a — !• , 
for a — > 



(4) 



a = is the Chern-Simons limit corresponding to the upper bound of the Gauss-Bonnet factor. When 

a — > 0, (H)) goes back to the Schwarzschild AdS soliton. For simplicity, in the following we will consider 
the Gauss-Bonnet factor in the range < a < I? jh. For the smoothness at the tip, we impose a period 
/? — f^^Zi-fr ^'^^ ^^'^ coordinate ip to remove the singularity. 

III. PHASE TRANSITION BETWEEN THE S-WAVE INSULATOR AND SUPERCONDUCTOR 



In the background of the Gauss-Bonnet-AdS soliton, we consider a Maxwell field and a charged complex 
scalar field coupled via the action 

1 



S = / d^xy/^ 



(5) 



Taking the ansatz of the matter fields as ^ — "tpir) and A — (j)(r)dt, we can get the equations of motion for 
the scalar field and gauge field (f) in the form 



/ r 



0-0, 



(6) 



(7) 



where the prime denotes the derivative with respect to r. 

In order to solve the above equations, we have to impose the boundary conditions at the tip r — Ts and at 
r — >■ oo. At the tip r = r^, the solutions behave as 

V' = V'o + i'lir - rs) + ^2{r - rsf H , 

^ = ^o+0i(r-r,) + 02(r-r,)2 + ... , (8) 

where V'i and (i = 0,1,2,- ••) are integration constants, and the Neumann-like boundary condition has 



been imposed to keep every physical quantity finite 



It is worth noticing that one can find a constant 



nonzero gauge field (j){rs) at r = rs, in contrary to that of the AdS black hole where (f>{r+) = at the horizon. 
Near the boundary r — > oo, we have asymptotic behaviors 

^=— + ^07''/' = /^-^' (9) 
where fi and p are interpreted as the chemical potential and charge density in the dual field theory respectively. 



Here A± = 2 ± •\/2 + m?L'^^. The coefficients il^- and "0+ both multiply normalizable modes of the scalar field 
equations and they correspond to the vacuum expectation values "0- O- >, '0+ 0+ > of operators 
dual to the scalar field according to the AdS/CFT correspondence. We can impose boundary conditions that 
either ')p- or 0+ vanish [5|, |6|. For simplicity, we will scale L = I and rg = 1 in the following just as in 

Before going further, we would like to give a comment. In the AdS soliton background, since at the tip 
(t>{rs) does not vanish in (6) and (7), which leads that the 0^ terms in the coupled equations cannot be got rid 
of as did in the AdS black hole case, we cannot count on the matching method to obtain the analytic result. 
Here we will apply the S-L method [23] to analytically investigate the properties of the s-wave holographic 
insulator/superconductor phase transition in the Gauss-Bonnet gravity. We will calculate the critical chemical 
potential to accommodate the phase transition and analytically derive the critical exponent of condensation 
operator. In addition, we will derive the relation between the charge density and the chemical potential near 
the phase transition point and examine the effect of the Gauss-Bonnet factor. 

A. Critical chemical potential 



Introducing a new variable z = rs/r, we can rewrite Eqs. (jG]) and ^ into 



where the prime here denotes the derivative with respective to z. 

It has been shown numerically that the solution is unstable and a hair can be developed when the chemical 
potential is bigger than a critical value, i.e., [i > ^c- For lower chemical potential, /i < ^c, the gravitational 
dual is an AdS soliton with a nonvanishing profile for the scalar field ip^ which can be viewed as an insulator 
phase 



18|. Thus, there is a phase transition between the insulator and superconductor phases around the 
critical chemical potential /Xc. 

At the critical chemical potential /ic, the scalar field -0 — 0. So near the critical point Eq. ((TT|) reduces to 



(12) 



With the Neumann-like boundary condition ([8|) for the gauge field at the tip r = r^, we can obtain the 
physical solution 0(z) — ji to Eq. (fT2|) when ji < fic- Considering the asymptotic behavior in Eq. ([9]), close 
to the critical point this solution indicates that p = near the AdS boundary 2 = 0, which gives fairly good 



agreement with numerical results in Ref. 



i8i. 



As ^ — ^ ^c, the scalar field equation (fTO|) reduces to 



f 1 



V' = 0. 



As in ^2S] . we introduce a trial function F{z) near the boundary z — which satisfies 

V'(z) ^ mz^^F{z), 



(13) 



(14) 



with i = + or i = — . Here we will impose the boundary condition F(0) = 1 and F'(0) = 0. Then, we can 
obtain the equation of motion for F{z) 

1 



F" + 



Defining a new function 



"2A, ^ 






■A.(A, - 1) 


z 






[ z^ 



z \ J zj z'^i 



F = 0. 



T{z) 



we can rewrite Eq. p5|) as 

{TF')' + T 



A,(A,-1) , A, //' 1 



z \ J z) z'^f 



F^O. 



(15) 



(16) 



(17) 



According to the Sturm-Liouville eigenvlaue problem j30[, we obtain the expression which can be used to 
estimate the minimum eigenvalue of /i^ 

,_J^'T{F'^-UF^)dz 



loVF'dz 



(18) 



with 

A,(A,;-1) , A, ff 1 



' z \ f z J z^f ' 

T 



In the following calculation, we will assume the trial function to be F{z) = 1 — az^, where a is a constant. 



In 



18| the condensate (0+) was numerically calculated in 5-dimensional Gauss-Bonnet AdS soliton back- 
ground simply by fixing V'- ~ 0. It has been shown numerically that the increase of the Gauss-Bonnet factor 
a results in the increase of the critical chemical potential, which means that the higher curvature correction 
will make it harder for the scalar hair to be condensated. Now we can use the S-L method to understand the 
condensation analytically. 

Using Eq. (1181) to compute the minimum eigenvalue of /i^ for i = -f , we can obtain the critical chemical 
potential /^^ for different strength of the curvature correction and the mass of the scalar field. As an example, 
we calculate the case of m^i^jj = in detail. From Eq. (flSl) . we obtain 

= S^^, (20) 



with 



^, , l-6a + \/l-4a(-l + 4a) 1 



3a3/2 16a2 



2V^(-3 + 20a) + 3(1 - Aaf log ) 

a/q;(1 - Aa 



3[-l + IQq - 30q2 + yi - 4a(l - 8a -f IQa^)] ^ 
(-10+15fl-6a2)7^ 

"(a, a) = ^ • (21) 

For different values of the Gauss-Bonnet factor, we can get the minimum eigenvalues of /^^ and the corre- 
sponding values of a, for example, ^^„j„ = 11.607 and a = 0.440 for a = 0.0001, /imj„ = 12.667 and a = 0.386 
for a ~ 0.1 and /ij„j„ ~ 14.365 and a = 0.257 for a ~ 0.2. Then, we have the critical chemical potential 
/^c = Mmin ,17]. In TableU we present the critical chemical potential fic for chosen values of the Gauss-Bonnet 
coupling a and various masses of the scalar field determined by fixing ra^L^f^. Comparing with numerical 
results, we find that the analytic results derived from S-L method are in good agreement with the numerical 
calculation. 

From TableU we observe that for the same mass of the scalar field, the critical chemical potential increases 
when the Gauss-Bonnet factor a becomes bigger. Our analytic result supports the observation obtained 

r 

numerically that the higher curvature correction can make the scalar hair more difficult to be developed pi 



22 



27 



31 



36{ . It is interesting to note that the analytical S-L method can give consistent critical chemical 
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TABLE I: The critical chemical potential /ic obtained by the analytical S-L method (left column) and from numerical 
calculation (right column) with chosen Gauss-Bonnet coupling and various masses of the scalar field for the s-wave 
holographic insulator and superconductor model. In order to compare with the results in Refs. [l^, [l^ , we also present 
the critical chemical potential for m^L^fj = —15/4. 



a 0.0001 0.1 0.2 



'2 T J- 


= 


3.407 


3.404 


3.559 


3.556 


3.790 


3.789 


2 r 2 


-1 


3.137 


3.135 


3.275 


3.272 


3.477 


3.475 


2 r 2 


-2 


2.817 


2.815 


2.937 


2.935 


3.106 


3.105 


2 r 2 


-3 


2.399 


2.396 


2.497 


2.494 


2.624 


2.622 


2 r 2 


15/4 


1.897 


1.888 


1.963 


1.960 


2.042 


2.039 



potential with the numerical result even when the mass of the scalar field is zero. This shows that the S-L 
method is more effective than the matching method. It was realized that the matching method cannot deal 
with the scalar field with zero mass 181, since for this case the Gauss-Bonnet term does not contribute to the 
analytic approximation in the matching method [181 119| . 

For the same strength of the curvature correction, with the increase of the mass of scalar field, the critical 



chemical potential fic becomes larger. This property also agrees well with the numerical result 18|. 

In the AdS black hole in Gauss-Bonnet gravity 26] , the S-L method was applied by choosing specific mass 

of the scalar field by fixing m^L^^j ~ —3. For choosing other nonzero values of the mass of scalar field, (4.16) 

there cannot be integrated analytically so that their (4.17), (4.18) and further steps cannot be derived. In 

the AdS soliton background, we found a much better situation. (18) above can be integrated analytically 

1 ^ 

when the mass of the scalar field satisfying the Breitenlohner-Freedman bound ^37,], so that analytically we 
can observe the condensation with the change of the mass of the scalar field. 

If we fix the scalar field mass by choosing values of m?L^ instead of m^L^g, the S-L method can give the 
same qualitative dependence of the critical chemical potential on the Gauss-Bonnet factor as described above 



when we study 



background 



the scalar operator (0+). This supports the numerical computation in the AdS black hole 



Now we concentrate on the scalar operator (O-) by imposing the condition -0+ = 0. Fixing the mass of the 
scalar field by choosing values of to^L^^, the scalar operator (O-) presents us qualitatively the same behavior 
of the condensation when we change the strength of the curvature correction as we observed above. With 
the increase of the Gauss-Bonnet coupling, the critical chemical potential will increase, which shows that the 
condensation will be harder to develop. But if we fix the mass of the scalar field by choosing values of m^L^, 
analytically we observed completely different condensation behavior as the Gauss-Bonnet coupling changes, 
see Table [ill Our analytical result got by using S-L method presents the same abnormal behavior as found 
numerically in 1^. Considering that choosing the mass of the scalar field by selecting the value of rn^L^g 
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contains directly the signature of Gauss-Bonnet factor in the scalar mass, we believe that this way of choosing 
the scalar field mass can disclose the correct consistent influence due to the Gauss-Bonnet coupling in various 
condensates. 

TABLE II: The critical chemical potential /ic obtained by using the analytical S-L method for the s-wave holographic 
insulator and superconductor model. The mass of scalar field is chosen by fixing m^L^ff ~ —15/4 and m^L^ = —15/4, 
respectively. 



a 


0.0001 


0.01 


0.05 


0.1 


m^L'iff = -15/4 


0.837 


0.839 


0.849 


0.862 


m^L'^ = -15/4 


0.836 


0.798 


0.648 


0.463 



B. Critical phenomena 

We will use the S-L method to analytically discuss the critical phenomena for the phase transition between 
the s-wave holographic insulator and superconductor in the Gauss-Bonnet gravity. We will concentrate on 
studying the critical exponent for condensation operator and the relations between the charge density and the 
chemical potential. 

The scalar field ip can be given by Eq. (|14p when —>■ /ic, so we can rewrite the equations of motion ([Tl]) as 

'/' 1\ 2(C',)2^2A.-4_p2 



J," 



= 0. (22) 



/ / 

Since the condensation for the scalar operator (Oi) is so small, we can expand 4'{z) in small (Oi) as 

^(z) ^ fic + {0^)x{z) + ■■■ . (23) 

Considering the boundary condition at the tip, we can get x(l) = and x'(l) = constant. After defining a 
function 



Jl + A(z^ - l)a - 1 
P{z) = ^ , 24 

we can obtain the equation of motion for xi^) 

(Px')' - m.)^^c J = 0. (25) 

According to the asymptotic behavior in Eq. Q , we can expand (j) when z — >■ as 

0(z) ^^i-pz^^^ic + m [x(0) + x'(o)z + ^x"(o)^' + •••]• (26) 

From the coefficients of the z° term, we can easily get 

fl-^l,^ (O.)x(O). (27) 
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If we set 

x{z) - 2{o,)^l,az), 

where the function ^(z) is the sohition to the following equation 



(28) 



/ 



(29) 



we will know that 



-yj(M-Mc)'/', 



(30) 



where ^(0) = ci — /o^[c2 + F{x)^x'^^^~^dx]jj^ with the integration constants ci and C2 determined by the 
boundary condition ^(z). For example, fixing m?L^^ = —15/4 and a — 0.0001, we can get ^(0) = 0.0815 



17|. 



when a = 0.330 which results in {0+) w 1.801(m - IJ^cf^- This agrees well with the result given in 

Note that our expression ([50]) is valid for all cases considered here, thus near the critical point, both of the 
scalar operators (0+) and (O-) satisfy {Oi) ^ [ji — /i^)^/^. This behavior holds for various values of Gauss- 
Bonnet couplings and masses of the scalar field. The analytic result supports the numerical computation 
[isi I27I that the phase transition between the s-wave holographic insulator and superconductor belongs to 
the second order and the critical exponent of the system takes the mean-field value 1/2. The Gauss-Bonnet 
coupling will not infiuence the result. 

Considering the coefficients of terms in Eq. p6|) . we find that x'(0) if z — 0, which is consistent 
with the following relation by integrating both sides of Eq. (|25|) 



VI - 4a - 1 Jo 



1 ,2Ai-4 



PF' 



f 



dz. 



(31) 



Comparing the coefficients of the term in Eq. 



we can express p as 



:(O.)x"(0). 



(32) 



From Eqs. ([25]) and (PT|) . we arrive at 
X"(0) 



X [z) 



P{z) 



y/l-Aa- 1 Jq 



1 ^2A,-4p_p2 



dz. 



(33) 



Using the above formula and Eq. (PU]) . we can deduce 



p = T{a,m){p. ~ p.c), 



(34) 
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where r{a,m) is a function of the Gauss-Bonnet couphng and the scalar field mass 

L {a,m) = = r-TT / 7 dz. (35) 



17|. 



Fixing m^L^fj = —15/4 and a — 0.0001, for example, we can get r{a,m) = 1.330 when a = 0.330, so that 
p — 1.330 (/I — fic) for considering the scalar operator (0+). This is consistent with the result given in 
Here we observed that the Gauss-Bonnet coupling will not alter the result. Our analytic finding of a linear 
relation between the charge density and the chemical potential p ^ (p — p.^) supports the numerical result 
reported in 1^ . 

IV. PHASE TRANSITION BETWEEN THE P-WAVE INSULATOR AND SUPERCONDUCTOR 

Since the S-L method is effective to obtain the s-wave holographic insulator and superconductor phase 
transition, we will use it to investigate analytically the p-wave holographic insulator and superconductor 
phase transition in the Gauss-Bonnet gravity which has not been constructed as far as we know. 

Considering an SU{2) Yang-Mills action in the bulk theory [38| 



^= / d^x^/^f-V F'^'^M, (36) 



4' 

where F^^, = d^A'^ - d^A^ + e'^'^'^A^^Al is the SU{2) Yang-Mills field strength and e"^" is the totally antisym- 
metric tensor with e^^'^ = -1-1. The A" are the components of the mixed- valued gauge fields A = A'Ji^T^'dx^, 
where t° are the three generators of the SU{2) algebra satisfy [t°,t''] = ^abc^c^ 

In order to construct a p-wave holographic insulator and superconductor in the Gauss-Bonnet gravity, we 



adopt the ansatz of the gauge fields as 
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39|-|41|, 



A{r) = (j>{r)T^dt + ^{r)T^dx. (37) 

Here we regard the U{\) symmetry generated by as the U{1) subgroup of SU{2). The gauge boson with 
nonzero component V'(^) along x-direction is charged under A^ = (/'(t'). According to AdS/CFT correspon- 
dence, (/)(r) and tpir) are dual to the chemical potential and the x-componcnt of some charged vector operator 
O in the boundary field theory respectively. The condensation of ipir) will spontaneously break the U{1) gauge 
symmetry and lead to a phase transition, which can be interpreted as a p-wave insulator and superconductor 
phase transition on the boundary. 

From the Yang-Mills action (l36|) . we can derive the following equations of motion 



7 + ^)*' 



*"+ l^ + r)*' + l7* = o> (38) 
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(39) 



where the prime denotes the derivative with respect to r. 

In order to solve the above equations of motion, we have to impose the boundary conditions for the (j){r) 
and tl){r) fields at the tip r = Vs and at r — oo. At the tip r = r^, the solutions have the same form just as 
Eq. (|H]) for the s-wave holographic insulator and superconductor model. But near the boundary r — )■ oo, we 
have different asymptotic behaviors 



V = Vo + , 



P 

P - —■• 



(40) 



where /i and p are interpreted as the chemical potential and charge density in the dual field theory, while V'o 
and 7/>2 may be identified as a source and the expectation value of the dual operator, respectively. Since we 
arc interested in the case where the dual operator is not sourced, we will set ipi) — Q and have a normalizable 
solution. 

A. Critical chemical potential 



Define the variable z = r^/r^ the equations of motion (j38|) and p9p can be expressed in the z coordinate as 



/ z 



(41) 



f z 



where the prime denotes the derivative with respective to z. 




(42) 




2.5 



3.0 



3.5 



4.0 



1.0 



1.5 



2.0 



2.5 



3.0 



FIG. 1: (color online) The condensates of the operator (O) = tp2 and charge density p with respect to the chemical 
potential fi for diflterent Gauss-Bonnet couplings a for the p-wave holographic insulator and superconductor model. 
The three lines from left to right correspond to increasing a, i.e., a = 0.0001 (red), 0.1 (blue and dashed) and 0.2 
(black) respectively. 
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Similar to the analysis in the previous section, if /i < fic, the field ip is nearly zero, i.e., tp 0. Thus, we 
can obtain the physical solution (l){z) — fi to Eq. (|42p when /i < fic. This is consistent with the numerical 
results in Fig. [1] which plot the condensates of the operator (O) — 'ip2 and charge density p with respect to 
the chemical potential /i for different Gauss-Bonnet couplings a. 

As p, ^ fic, Eq. (liTI) will become 



(43) 



We can also define a trial function F{z) near the boundary z ~ just as in the last section 



i;iz) ~ {0)z^F{z), 



(44) 



with the boundary condition F{0) ~ 1 and -F"(0) ~ 0. Therefore, the equation of motion for F{z) is given by 



F" 



Introducing a new function 



z f 



T{z) 



zf z^f 



(45) 



+ 4(z4 - l)a - 1) 



(46) 



we can rewrite Eq. p5)) as 



iTF'y+T[±+^J^ + -^]F^O. 



(47) 



Defining the following parameters 



4 2/' T 



we find that, following the Sturm-Liouville eigenvlaue problem 
obtained from variation of the following functional 
.^IoT{F''-UF^)dz 



(48) 



the minimum eigenvalue of can be 



J^VF^dz 

12q;(1 — 4a) + a(8 — 3a + 20aa) - 8a(l - Aaf^'^ 



4[6 + (3a- 8)a]aS/2 
+3(1 - 4Q!)[4a + a^(4a - 1)] log 



(1 - Aa)a 



(49) 



_(2a + yS)2_ 

where we have assumed the trial function to be F{z) = 1 — az^ with a constant a. 

We can easily obtain the minimum eigenvalues of and the corresponding values of a for different Gauss- 
Bonnet couplings a, for example, /imi„ = 5.140 and a = 0.338 for a = 0.0001, ~ 5.697 and a ~ 0.305 for 
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a = 0.1 and — 6.691 and a = 0.223 for a = 0.2. Thus, we get the critical chemical potential Hc — /i„ 



17 1 which has been shown in Table Hill for fixed value of the Gauss-Bonnet factor a. In order to compare 
with numerical results, we also give the critical chemical potential obtained by using the shooting method. 
Obviously, the agreement of the analytic results derived from S-L method with the numerical calculation is 
quite impressive. 

TABLE III: The critical chemical potential Hc obtained by the analytical S-L method (left column) and from numerical 
calculation (right column) with fixed Gauss-Bonnet coupling for the p-wave holographic insulator and superconductor 
model. Note that our result reduces to the result in Ref. [l6] if q — > 0. 



a 


0.0001 


0.1 


0.2 


He 


2.267 2.265 


2.387 2.385 


2.587 2.585 



From Table IIIIl we also find that the critical chemical potential increases when the Gauss-Bonnet factor a 
becomes bigger, which shows that the higher order curvature corrections in general make the condensation 
harder to form, just as observed for the s-wave holographic insulator and superconductor model. This property 
agrees well with the numerical result shown in Fig. [1] 



B. Critical phenomena 

With Eq. (|44l) . when fi — > /ic the equation of motion ()42|) can be rewrited as 

1\ {Ofz^F^ 



.... , .=0. (50) 

f z) f 

Note that the condensation value of "ipiz) is so small, we will expand (l){z) in small (O) as 

</)(z) ~ Me + (0)xW + • • • , (51) 

with the boundary condition = at the tip. Using the function defined in Eq. ([M]). we can get the 
equation of motion for x(z) 

{Px'Y ~ {0)fic^-^ = 0. (52) 

Near z — 0, we can also expand (f> as 

0(z) :^fi-pz^^ fi, + {O) [x(0) + x'(0)z + ^x"(0)z' + •••]. (53) 
Comparing the coefficients of the term in both sides of the above formula, we can obtain 

- ~ (O)X(O). (54) 
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Considering the following equation for ^(z) 



with 



X{z) = {0)fi,^{z), 



(55) 



(56) 



we will have 



{0) = 



tt7j(m-Mc) 



1/2 



(57) 



wmv 

where ^(0) = ci — Jq[c2 + F{x)'^x^dx]jj^ with the integration constants ci and C2 determined by the 
boundary condition x{z). For example, for the case of a = 0.0001, we can get ^(0) = 0.0673 when a = 0.338 
which results in (O) « 2.560{fi - fic)^/^ . This is in good agreement with the result given in |16l . Il7| . 



It should be noted that the relation ([57)1 is valid for all cases considered here, so the condensation (O) ^ 
(/i — /ic)^/^ near the critical point for various values of Gauss-Bonnet couplings, which agrees well the numerical 
results in Fig. [T] that the phase transition between the p-wave holographic insulator and superconductor 
belongs to the second order and the critical exponent of the system takes the mean-field value 1/2. 

From the coefficients of the term in Eq. ((53|) . we obtain that x'(0) — > which is consistent with the 
following relation by making integration of both sides of Eq. ([52|) 



-dz. 



For the coefficients of the term in Eq. ([5^ . we have 

p = -^(O)x"(0) = r(a)(M-M,), 

where T{a) is only the function of the Gauss-Bonnet couplings which can be given by 



(58) 



(59) 



r(a) = 



-dz. 



(71^4^- 1)^(0) 7o / 
For example, we can obtain T{a) ~ 1.126 when a = 0.338 for a = 0.0001, i.e., the linear relation p 



(60) 



17[. Still we notice that the Gauss-Bonnet coupling 



1. 1 26 (/i — /ic), which agrees well with the result given in 
will not change the result. The analytic finding of a linear relation between the charge density and the chemical 
potential p ^ (p — pc) is consistent with the numerical result presented in Fig. [1] 
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V. CONCLUSIONS 

We have applied the S-L method to investigate analytically the condensation and critical phenomena of 
the phase transition between the holographic insulator and superconductor in the Gauss-Bonnet gravity. We 
found that unlike the analytic matching method, the S-L method is effective to obtain the analytic rcsiilts in 
the AdS soliton background both for s-wave (the scalar field) and p-wave (the vector field) models. Different 
from the AdS black hole in the Gauss-Bonnet gravity, in the AdS soliton spacetime we observed that the 
S-L method can bring us results of condensation for different values of the scalar field mass satisfying the 
Breitenlohncr-Frccdman bound. For the massless scalar field, the information of the Gauss-Bonnet coupling 
can still be kept in the S-L method. With this analytic method, we also found that it is more appropriate 
to choose the mass of the scalar field by selecting the value of m^Lgjf . The analytic results derived from the 
S-L method for the s-wave and p-wave holographic insulator /superconductor phase transitions support the 
numerical computations and show that the higher curvature corrections make it harder for the condensation 
to form. 

Furthermore, comparing with the matching method, we found that in the AdS soliton in the Gauss-Bonnet 
gravity the S-L method can present us analytic results on critical exponent of condensation operator and the 
relation between the charge density and the chemical potential near the phase transition point. We observed 
that effect of the Gauss-Bonnet factor cannot modify the critical phenomena. The analytic results can be 
used to back up the numerical findings in both s-wave and p-wave insulator /superconductor models of the 
Gauss-Bonnet gravity. 
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